Bose-Einstein condensation into non-equilibrium states studied by condensate focusing 
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We report the formation of Bose-Einstein condensates into non-equilibrium states. Our conden- 
sates are much longer than equilibrium condensates with the same number of atoms, show strong 
phase fluctuations and have a dynamical evolution similar to that of quadrupole shape oscillations of 
regular condensates. The condensates emerge in elongated traps as the result of local thermalization 
when the nucleation time is short compared to the axial oscillation time. We introduce condensate 
focusing as a powerful method to extract the phase-coherence length of Bose-Einstein condensates. 
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Among the quantum fluids, the quantum gases are es- 
pecially suited to study the kinetics of Bose-Einstein con- 
densation (BEG). The modest density of the quantum 
gases makes it possible to stretch the time of condensate 
formation to values allowing detailed experimental inves- 
tigation. In external potentials condensates have a char- 
acteristic equilibrium shape, known as the Thomas-Fermi 
shape, that differs markedly from the shape of thermal 
clouds. Since the pioneering experiments on BEG in the 
alkali systems H, |^, this feature has been observed 
by many groups and is used routinely to measure the 
condensate fraction jj]. Phase coherence is another key 
property of equilibrium condensates and was established 
in the first interference experiments [|[ |[ 0, ||. 

Equilibrium condensates ^ are produced by quasi- 
static growth, where heat extraction limits the formation 
rate. The condensate nucleates as a small feature in the 
center of the trap and grows as long as heat is extracted 
from the sample. To observe the formation kinetics, the 
gas has to be brought out of equilibrium, in practice by 
shock cooling. Since the first experiment on condensate 
growth, by Miesner et al. |^ , this is done by fast RF re- 
moval of the most energetic atoms from the trap. Start- 
ing from a thermal gas just above the phase transition 
temperature {Tc), the condensate appears as the result 
of thermalization. Miesner et al. |^ observed the growth 
under adiabatic conditions. Kohl et al. |0 continued the 
extraction of heat and atoms, also during growth. In 
both experiments, the condensate was observed to grow 
from the center of the trap, like in the quasi-static limit. 

Kagan et al. [ pT| pointed out that qualitatively differ- 
ent stages have to be distinguished in the formation of 
equilibrium condensates with a large number of atoms. 
The early stage (kinetic stage) is governed by Boltzmann 
kinetic processes and leads to a preferential occupation 
of the lowest energy levels. Once a substantial fraction 
of the atoms gathers within an energy band of the or- 
der of the chemical potential of the emerging conden- 
sate during formation, their density fluctuations are sup- 
pressed in a fast interaction-dominated regime governed 



by a non-linear equation for the boson held. The appear- 
ing phase-fluctuating condensate then grows and the con- 
densed fraction approaches its equilibrium value. How- 
ever, the phase fluctuations still persist, giving rise to 
dynamically evolving flow patterns in search for the true 
equilibrium state. In elongated 3D trapped gases the 
phase fluctuations can be pronounced even under equi- 
librium conditions as was predicted by Petrov et al. 
and observed experimentally by Dettmer et al. [ p^ . 

In this Letter we report the formation of condensates 
into non-equilibrium states and a new path towards equi- 
librium in elongated traps. In contrast to the previ- 
ous experiments our results were obtained starting from 
thermal clouds deep in the cross-over regime to hydro- 
dynamic behavior. The condensates are much longer 
than equilibrium condensates with the same number of 
atoms. Moreover, they display strong phase fluctua- 
tions and a dynamical evolution similar to that of a 
quadrupole shape oscillation decaying towards equilib- 
rium. We identify l/uJz as a characteristic time that 
should be addressed explicitly for elongated cylindrical 
harmonic traps, i.e. for traps with ujp^ ujz, where ujp and 
LOz are the radial and axial angular frequencies, respec- 
tively. We show that these exotic condensates emerge 
as the result of local thermalization when the nucleation 
time is short as compared to l/a;^. The dynamical evo- 
lution of the condensate in the trap has to be dealt with 
explicitly to properly interpret time-of-flight absorption 
images. In this context we introduce condensate focus- 
ing as an alternative to Bragg scattering 1^ for measur- 
ing the phase-coherence length of phase-fluctuating Bose- 
Einstein condensates. 

In the previous experiments on condensate formation 
the phase fluctuations were not studied. The results of 
Miesner et al. were compared to an analytical expres- 
sion for adiabatic growth of a condensate from a thermal 
cloud, derived by Gardiner et al. |Q . Although a qualita- 
tive agreement between theory and experiment was read- 
ily obtained it turned out to be impossible to obtain de- 
tailed agreement at the quantitative level |n5, 16 . In the 
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experiment of Kohl et al. [ |lO| quantitative agreement with 
the quantum kinetic approach (see refs. |l5|, ^) was 
obtained for strong truncation, whereas for weak trunca- 
tion the observed behavior differed distinctly from theory. 

In our experiments we typically load 4 x 10^ atoms of 
^^Wo in the |5^S'i/2, ^ = 2, mp — 2) state into a horizon- 
tal loffe-Pritchard quadrupole trap with ojp = 27tx 477(2) 
Hz and = 27r x 20.8(1) Hz (no RF dressing) The 
trap minimum Bq = 88.6(1) /iT corresponds to a radio 
frequency of i^q — 620 kHz as calibrated against atom 
laser output coupling jl8|. The trap minimum shows a 
long term drift of 5 kHz/hr. At full current thermal drift 
effects are less than 1 kHz/s. 

The gas is prepared by forced RF evaporation at a 
final rate of !>tr = —433 kHz/s to a value I'tr.a = 740 kHz, 
followed by 20 ms of plain evaporation at J^tr.a- As 

|i>tr/(t'tr,a ~ t^o)| < ^^z (1) 

this yields a static, purely thermal cloud of A^i « 5 x 10^ 
atoms at a temperature Tq = 1-3(1) /iK as determined 
from the axial size Iz = [2kTo/m,u)'^]^^'^ of the cloud in 
the trap, measured by time of flight absorption imaging 
shortly {t < 1/uJz) after release from the trap. We calcu- 
late a central density tiq ~ 4 x 10^** cm^'^, corresponding 
to a mean free path Ao = (2^/^rtoCT)^^ ~ 3 fim and a 
collision rate t^~J « noVth<^ « 5000 s^^ [|l9|. Radially 
we find \o/lp ~ 0.5. Axially we have Xq/Iz ~ 0.02 to be 
compared with the values Xo/lz ~ 0.5 and 0.3 used in 
previous experiments ||] and ||l^ on condensate forma- 
tion. Hence, our thermal samples are prepared far deeper 
in the cross-over to the hydrodynamic regime. 

The hydrodynamic behavior manifests itself in the 
damping time tq and a shift of the frequency ujq of the 
quadrupole shape oscillation as well as in an anisotropic 
expansion of the cloud after release from the trap. We 
measured uJzTq = 10(3) and uq/uJz — 1.56(5). For 
i^zTQ = 10(3) theory predicts luq/lUz = 1.57(2) | |2C| ] 
and ujq/uJz = y/l2/5 w 1.55 for the hydrodynamic 
limit in very elongated traps |2^, Further, writ- 

ing f3 = ujz/i^p, from the scaling theory ||2^, we find 
in this limit Vz/vp w 2.6/3 w 0.11 for the ratio of axial 
to radial expansion velocity after release from the trap. 
We measured Vz/vp = 0.78(2) and a final temperature 
Too = 0.94(4) /iK as determined from the axial expan- 
sion, which implies that the gas expands hydrodynami- 
cally and cools only briefly before the expansion becomes 
ballistic. For isentropic hydrodynamic expansions the de- 
generacy parameter is conserved: n{T)/no = {Too/ToY^^- 
Then, assuming that the expansion is hydrodynamic ini- 
tially and becomes ballistic after a time Tfreeze, the scaling 
theory (2ll gives a relation between Tq/T^c and Vz/vp. In 
elongated traps, for Vz/vp ^ (3, we find 

To/Too « [1 + 2{vp/vzf]/i. (2) 

The observed ratio Vz/vp and temperature Tq lead to 
Too = 0.91(7), which coincides with Too measured 



from the axial expansion. From the obtained value of 
n{Ti-ccc7.o) / no we calculate rtrcozo ~ 0.3 ms, i.e. a cou- 
ple of collision times. This is consistent with the ra- 
tio Xo/lp ~ 0.5 mentioned above and confirms the pic- 
ture that we are operating near the onset of hydrody- 
namic behavior in the radial direction. For complete- 
ness we verified that the expansion becomes isotropic, 
Vz/vp — 1-02(4), when the atom number is reduced by a 
factor of 30. 

Once the thermal cloud is prepared we distinguish 
three distinct stages. First, in the truncation stage, the 
radio frequency is set to the value i/tr,b = 660 kHz. This 
stage has a duration = 1 ms, which is chosen to be long 
enough (^tr > 1/wp) to allow atoms with radial energy Sp 
larger than the RF truncation energy etr to escape from 
the trap, yet is short enough to disallow evaporative cool- 
ing. We found that in this stage 50% of the atoms are 
removed. Notice that due to the finite radial escape hori- 
zon {\o/lp « 0.4) the ejection is not expected to be com- 
plete. Futhermore, the escape efficiency is anisotropic as 
a result of gravitational sag. The truncation energy etr 
covers the range 3 - 5 /xK depending on the position of 
the truncation edge in the gravity field and is lowered by 
an additional 1 /xK due to RF-dressing (Rabi frequency 
fiif w 27r X 14 kHz). At the start of the second stage, 
the thermalization stage, the radio frequency is stepped 
back up for a time tth to the frequency iytr,& to allow the 
gas to thermalize under formation of a condensate. The 
value I'tr.a is chosen to eliminate any appreciable evap- 
orative cooling. The third stage, the expansion stage, 
starts by switching off the trap and covers the time of 
flight T after which the sample is absorption imaged on 
the \5^Si/2,F = 2) ^ \5^P3/2,F = 1, 2 or 3) transition. 

To follow the evolution of the trapped gas after the 
truncation we took time-of-flight absorption images for 
a range of evolution times t = ttr -I- tth and a fixed ex- 
pansion time r. The images show a bimodal distribution, 
indicating that the truncation procedure results in BEG. 
The condensate fraction grows to a final value of 6% with 
a characteristic time of 6 ms. This corresponds to SOtcoi, 
in accordance with previous experiments. 

Rather than discussing the details of the growth ki- 
netics we emphasize that our condensates nucleate into 
non-equilibrium states. In Fig.l we plot the Thomas- 
Fermi half-length Lz obtained with the standard fitting 
procedure of a bimodal distribution to our data Q. For 
the shortest expansion time, r = 2.8 ms, the axial size 
of the condensate image equals to good approximation 
(t <C 1/wz) the axial size of the condensate in the 
trap. We see that Lz{t) is initially oversized by a fac- 
tor Lz{0)/ Lz{oo) = 2.2(3) and rapidly decreases to reach 
its equilibrium size after roughly one strongly damped 
shape oscillation (see open triangles Fig.l). Hence, the 
condensate is clearly not in equilibrium p3| . 

The formation of oversized condensates follows from 
the local formation concept underlying ref . pl| . Starting 
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FIG. 1: Condensate length Lz{t,T) versus evolution time t 
for three different expansion times r. The dark grey line is 
a guide to the eye. The light grey line represents a fit to a 
damped quadrupole shape oscillation. 



with a thermal cloud of A'i = 5 x 10*^ atoms in an 11.5 
mK deep harmonic trap at a temperature Tq = 1.3 fiK, 
we calculate that 55% of the atoms remain trapped after 
all atoms with energy e > Stv = 3.4 /iK are removed. The 
gas will rethermalize well within a time short compared 
to I/lOz- Hence, the resulting temperature varies along 
the trap axis. In this respect the thermalization is a local 



phenomenon. The local Tc is given by |24 



kTc{z) K. 1.28 liiOp [ni-o{z)r^ 
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(3) 



where Vp = [H/mujpY^^ is the radial oscillator length and 
niuiz) the atom number per unit length at position z. 
We find that the local temperature T{z) is lower than the 
local Tc{z) over a length of order Iz- In view of the sim- 
plicity of this model we consider this as good qualitative 
agreement with experiment. 

To further investigate the formation process we intro- 
duce condensate focusing. In our case one-dimensional 
focusing results from axial contraction of the expanding 
cloud when the gas is released from the trap during the 
inward phase of a shape oscillation. The focus is best 
demonstrated by plotting the axial size Lz{t,T) of the 
condensate as a function of expansion time r after a fixed 
evolution of t = 11 ms in the trap (see Fig. 2). The axial 
size is seen first to decrease and to increase again later 
as expected for a focus [|6) . 

We elucidate the focusing in relation to the low- 
frequency branch of the m = quadrupole shape os- 
cillation of the Thomas-Fermi condensate. In the limit 
of linear response the scaling parameter bz for the ax- 
ial size of the condensate in the trap can be written as 
bz{t) = Lz{t)/Lz{0) = 1 -|- Qz cos ujQt, with Uz <^ 1 the 
scaled axial amplitude and t the evolution time. Releas- 
ing the gas at time t and observing the cloud after an 
expansion time r 3> l/^^p, the scaling theory |21, |2^ 



offers an approximate expression for the scaled axial size 



bz{t,T) sa 1 -f {tt(3u!z/2 - a^WQsinwgt)? 



(4) 



FIG. 2: Condensate length Lz{t,T) versus expansion time r 
after 11 ms of evolution (10 point averages). The grey 
(dashed) line is a fit to the scaling equations including (ex- 
cluding) phase fluctuation broadening. Inset: axial density 
profile showing phase fluctuation stripes ai t — 6.8 ms. 



The first term under the brackets corresponds to the axial 
expansion kick caused by the declining chemical potential 
at trap release. The second term represents the scaled 
axial dilatation velocity in the trap at the moment of 
release. From Eq.(H) we see that the axial dilatation 
field reaches a real (virtual) focus for positive (negative) 
values of the expansion time Tfocus = (azWgsinwgt — 
^(3uJz)~^- A real focus can be obtained already for small 
amplitudes, Gz > 13. 

For r ^ the radial expansion is described by 

bp{t,T) « iLJpT{l — jGz cosujqt), i.e. shows no focus. As 
the radial size remains finite, the chemical potential will 
build up near the focus until the compression is balanced 
and the axial size starts to expand again. As follows from 
the scaling equations (see ref. |2^ or [0), at maximum 
compression the axial size is reduced by a factor of 

Our data match a scaled focal size bz{t, Tfocus) = 
0.49(6), i.e. the focus is strongly broadened as compared 
to the minimum scaled size The broadening is at- 
tributed to local variations in expansion velocity caused 
by the presence of phase fiuctuations in our condensates. 
After some expansion these variations give rise to irregu- 
lar stripes (see the inset of Fig. 2) as previously observed 
in Hannover ||l^. At the focus, the axial distribution 
maps linearly onto the momentum distribution due to the 
phase fiuctuations in the original condensate. The scaled 
focal size is given by 6z(i, Tfocus) {fi/mL^)Tfocus/Lz, 
where is the phase coherence length and {h/mL^) 
characterizes the expansion velocity due to the phase 
fluctuations. For equilibrium phase fluctuations close 
to Tc we estimate Lz/L^ — 7(4) (see ref. Q) and 
bz ~ 0.05. As the observed focal size is larger by an 
order of magnitude, the phase fluctuations have a non- 
equilibrium origin with a phase coherence length of only 
^ I fim. The decay of the phase fluctuations as a 
function of evolution time is subject of further investiga- 
tion in our group. 
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We point out that for a thermal cloud driven on 
the low-frequency quadrupole mode the phenomenon of 
focusing is absent except deeply in the hydrodynamic 
regime. For a small ratio T^o/Tq the scaled focal size 
will be given by ^^(t, rfocus) ~ Vth{Toc)T{ocus/lz{To) « 
(roo/To)i/2c^,rfocus < 1. 

In Fig. 1 we also show the oscillation in the axial size of 
the condensate as observed for 15.3 ms (open squares) 
and 25.3 ms (crosses) of expansion Due to en- 

hancement by focusing the amplitude of the oscillation 
has increased as compared to the 2.8 ms results. For 
T = 25.3 ms the shape oscillation is seen to exceed 
the noise for at least 100 ms. This oscillation can be 
described by a linear response expression for evolution 
times t >20 ms where < 0.2. We measure a damping 
time of TQ = 50(9) ms and a frequency ratio ujq/uJz — 
1.54(4). The latter is slightly lower than the frequency 
expected for a quadrupole shape oscillation of a pure con- 
densate in very elongated traps, coq/lOz ~ ^5/2 « 1.58 
pl| . A 5% negative frequency shift was observed for 
the quadrupole mode in Na condensates just below Tc 
pof and is consistent with theory |30[ |. 

The condensate focusing is a particularly useful tool. It 
separates the condensate from the thermal cloud in time- 
of-flight absorption imaging. This enables the obser- 
vation of small condensate fractions and measurements 
close to Tc- It also allows the observation of small ther- 
mal clouds and thus extending of the dynamic range of 
time-of-flight thermometry. The focal size and shape al- 
low to quantify the phase fluctuations inside condensates. 
Extended to two dimensions the focusing can serve imag- 
ing applications, with the mean field of the quantum fluid 
acting as a tunable component in atom optics. 
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